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Abstract 

We consider the long term fate and evolution of cold degenerate 
stars under the action of gravity alone. Although such stars cannot 
emit radiation through the Hawking mechanism, the wave function of 
the star will contain a small admixture of black hole states. These 
black hole states will emit radiation and hence the star can lose its 
mass energy in the long term. We discuss the allowed range of possible 
degenerate stellar evolution within this framework. 
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I. INTRODUCTION 


White dwarfs are usually considered to be the dual state of stellar evolu¬ 
tion for nearly all stars in the universe [1]. The remaining stars are destined 
to end their lives as neutron stars or black holes. A fundamental astrophysi- 
cal question is to determine the long term fate and evolution of these stellar 
remnants. If the proton is unstable through some process at the unification 
scale Mx [2], then long term stellar evolution will be governed by the physics 
of proton decay with a time scale of 

Tp ~ 10^V(Mx/10^®GeV)^ . (1) 

This evolutionary path has been studied previously [1, 3, 4]. In this paper, 
we consider the ultimate fate of these cold degenerate stars in the absence 
of unihcation scale proton decay, i.e., due to the influence of only the gravi¬ 
tational force. In particular, we consider processes involving black holes and 
their effects on the long term evolution of these degenerate stars. 

If the proton were truly stable, then nature presents us with a curious 
state of affairs. Black holes of stellar mass, which are much more tightly 
bound than degenerate stars, will evaporate through the Hawking effect [5] 
with a lifetime of only ~ 10®® yr [6-8]. Although this time scale seems long 
compared to the current age of the universe (~ 10^® yr), it is vastly shorter 
than the expected lifetimes of white dwarfs and neutron stars. These stellar 
remnants are essentially zero temperature objects; they exist in hydrostatic 
equilibrium due to quantum mechanical degeneracy pressure [9]. In the ab¬ 
sence of proton decay, these remnants would live almost forever [1, 10]. In 
this case, black holes would evaporate almost instantaneously compared to 
white dwarfs and neutron stars. 

Since degenerate stars are far too extended to have event horizons, they 
cannot emit Hawking radiation in the usual manner available to black holes 
[11] (see also [6, 8, 12, 13]). Ref. [11] proved that in a globally static 
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spacetime, such as that produced by a degenerate star, there is no particle 
creation. These results break down in black hole spacetimes because they are 
not globally static, where static means both time independent and invariant 
under time reversal. For black holes, the staticity breaks down at and beyond 
the event horizon. 

After cooling to essentially zero temperature, degenerate stars can only 
lose energy through processes that utilize black holes, which can radiate. This 
avenue is studied in this present paper. Notice that this approach implicitly 
includes energy loss from proton decay via gravitation [14-17], since such 
processes involve microscopic black holes. 

Motivated by the above discussion, we write the wave function of a cold 
degenerate star in the suggestive form 

|star) = cos^liV*) + sin6\BH ), (2) 

where |iV*) represents the usual stellar conhgurations (superpositions of quan¬ 
tum mechanical states of TV* = few xlO®^ particles) and \BH) represents 
superpositions of the possible black hole states. The total probability of the 
star being in a black hole state or containing a black hole is thus V = sin^ 9 ~ 

9 ^ < 1 . 

In general, the black hole part of the wave function will contain many 
different contributions, i.e.. 


|i?//) = ^A,|Mbh), (3) 
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where |Mbh) represents the wave function for a black hole of mass Mbh and 
must be a superposition of many states. The quantity \ Aj\‘^ is the correspond¬ 
ing probability of the star being in a black hole state. The total probability 
is thus given by 

V = Ar? e = ■ W 
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We can estimate the probabilities Pj = from several different mecha¬ 
nisms, as discussed below. 

II. BLACK HOLE PROCESSES IN STARS 


A. Dyson Tunneling 

On sufficiently long time scales, white dwarfs and neutron stars can ex¬ 
perience tunneling events in which the star, or portions of the star, tunnel 
through the potential energy barrier produced by degeneracy pressure [10]. 
This process can, in principle, result in the formation of black holes. 

We hrst consider the probability of tunneling through the barrier provided 
by degeneracy pressure. Such tunneling is a necessary but not sufficient 
condition for black hole formation. The time scales for tunneling processes 
are given by the standard formula. 


r = To e^"^ , (5) 

where tq is the natural oscillation time of the system. The action integral St 
can be written 

St = ^ J[2MV{r)f^dr, ( 6 ) 

where V (r) represents the height of the potential energy barrier as a function 
of the radial coordinate r. 

For a white dwarf supported by degenerate electrons, the potential energy 
can be obtained from integrals over the distribution functions for fermions [7]. 
Following standard formalism, we dehne x = pp/me, where pp is the Fermi 
momentum and rrie is the electron mass. The number density of electrons is 
then given by Ue = x^ml/3n'^ (we take h = 1 = c). We want to consider a 
small spherical volume, of radius r, within the star that contains N electrons 
and {A/Z)N nucleons. After some algebra, the potential energy Vd[r(a;)] due 
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to degeneracy pressure can be written in the form 


Vd = NrUe-^-^ {1 + 2x'^) — ln[a;+ (1 + — 1 


= Nrue^ix) , 

( 7 ) 

where the second equality dehnes a dimensionless function <f)(a:). In the limit 
of small X (low density), this potential has the usual form associated with 
nonrelativistic degenerate fermions, i.e., $ oc oc r~^. In the opposite (rel¬ 
ativistic) limit of large x, <h oc a: oc r~^. We also must include an additional 
contribution to the potential energy due to gravity. This potential can be 
represented with the approximate form 




-/ 


r — rs ^ 


( 8 ) 


where = ANmp is the total mass of the volume, rg = 2 GMn is the 
corresponding Schwarzschild radius, and we have dehned a geometrical factor 
/ ~ 3/5. Unless the radius r approaches the Schwarzschild radius rg, the 
degeneracy part of the potential completely dominates that of gravity. As 
a result, to a good approximation, we only need to include Vd in the action 
integral 0, i.e.. 


ST^2{QTi/Af/^[2{A/Z)mp/meY/^N^'^ / —^^/^(x). (9) 

J xq 

The integration ranges from starting values xq near unity, corresponding to 
the initial density of the volume, to a large upper cutoff hi ^ 1, corresponding 
to black-hole-like densities. The exact value of the upper cutoff does not 
matter in practice because the integral has almost all of its support for small 
X (low density). 

For white dwarfs, with A/Z ^ 2, we have 


St M335Af‘‘''V(io), 


( 10 ) 


where /(xq) is the integral evaluated for a given starting value of Xq, i.e., a 
given starting density. The fastest tunneling rate will occur in the center of 
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the star where the density (and xq) is maximal. For a typical white dwarf of 
mass 0.5 Mq, the central density is abont 10® g/cm^, xq ~ 0.80, and /(xq) ~ 
1.56. For a more massive white dwarf of 1.0 Mq, the central density is abont 
3 X 10'^ g/cm^, xo ~ 2.5, and /(xq) ~ 1.0. 

For any given nnmber N of particles that are reqnired to form a black 
hole, eqnation determines the tnnneling probability oc exp[—S' t]. We 
can thns estimate the probability (per natnral oscillation time tq) for white 
dwarfs to prodnce black holes throngh a tnnneling process. 


p = g-335/(xo)7V4/3 


( 11 ) 


where N determines the mass of the black hole (Mbh ~ {A/Z)Nmp). This 
process is exponentially snppressed as the nnmber N of particles increases, 
so the process is dominated by the minimnm valne of N. 

As a reference point, if the minimnm mass for a black hole is given by the 
Planck scale, Mpi ~ 2 x 10“® g, we obtain N ~ 10^® and hence Sp ~ 10^^. 
For this incredibly large valne for the action, the natnral vibration time scale 
To becomes completely irrelevant. 

The fastest possible rate for this process would be for the case N = 2, 
and the action would be Sp ~ 844 I (xq) . In the center of a relatively massive 
white dwarf, where / ~ 1, the corresponding time scale is r 
10336 yp This time scale is still long compared to that of gravitationally 
induced proton decay [1, 14-16, 18]. We thus conclude that the tunneling 
of aggregates of baryons into black holes will not be an important process 
within white dwarfs. 

The result is similar for the case of neutron stars [10]. In this case, 
we must replace the electron mass rrie with the nucleon mass in all of the 
formulae and use A/Z = 1. The resulting action is much smaller, 

St ~ 5.4iV^/3/(xo). (12) 

At the center of a neutron star, the density is ~ 10^® g/cm^ and hence 
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approaches nuclear denisty. As a result, xq ~ 1, and hence /(xq) ~ 1.4. 
The action then becomes St ~ For the fastest case of only two 

nucleons, N = 2, we have St ~ 19. This small action would seemingly imply 
a reasonably high probability of tunneling into a black hole. However, this 
calculation only takes into account the tunneling of neutrons through the 
potential barrier provided by their degeneracy pressure. For small values of 
N and high densities, the neutrons start to overlap, and the quark structure of 
the nucleons must also be taken into account. The process of small numbers 
of quarks tunneling into a black hole is one allowed channel for gravitationally 
induced proton decay and the estimated time scale is about 10^^ yr (see [14, 
18] and also the following discussion). 

B. Kinematic Stellar Model 

In this approximation, we consider the star to be a collection of particles 
of mass m with a number density n. We want to calculate the probability 
that N particles happen to lie within their own Schwarzschild radius, i.e., 
that of a black hole of mass Mbh = mN. In natural units, the volume of the 
hypothetical black hole is given by 

I4h = (13) 

These N particles generally occupy a much larger volume Vq given by 

Vo = N/n. (14) 

Thus, the probability pi that a single particle lies within the black hole 
volume is simply 

Pi = ^ = , (15) 

ho 

and the probability pn that N particles he within the black hole volume is 

Pn = {pi)^ = {8nm^N^Mpi~^)^. (16) 
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Since the star contains N^/N volumes of N particles, the total probability 
Pn that the star contains a black hole composed of N particles is given by 

Pn = N^N-\8nm^N^Mpi-^)^ = . (17) 

In the second equality, we have dehned A = Snm^MpP^ ^ 1, 

where the numerical value assumes a nuclear density to evaluate n. For a 
white dwarf density, n will be far smaller, and the value of A will decrease 
accordingly. For all values of N less than iV*, the probability is a decreasing 
function of N. In other words, the most probable black hole to be formed 
kinematically will have the smallest possible value of N. 

For the case in which the minimum black hole mass is the Planck mass, 
N ^ 6 X 10^®, and hence — ln[P/v] ~ 10^^. Loosely speaking, this time 
scale is somewhat longer than, but comparable to, the tunneling time scale 
calculated previously. On the other hand, for small values of N, the time 
scale is Iqsq Pqj, g] 2 iall values of N, one can also construct 

a kinematic model of the proton as a system of iV = 3 quarks and make a 
similar calculation to model proton decay [18]. 

C. Virtual Black Holes 

In any quantum theory, one expects to hud vacuum fluctuations asso¬ 
ciated with the fundamental excitations of the theory. Thus, in electrody¬ 
namics, one has the possibility of forming electron-positron pairs for a short 
time directly out of the vacuum. The existence of such processes can be 
observed indirectly by many quantum phenomena, such as the Casmir effect. 
However, this example illustrates only the simplest such possibility. One 
must also include in the vacuum processes all excitations of the theory. One 
should thus include the possibility of the production of proton-antiproton 
pairs, or even monopole-antimonopole pairs. These processes will generally 
be highly suppressed relative to the electron-positron amplitudes by virtue 
of their correspondingly large masses. 



In gravitation, one therefore expects not only to find virtual gravitons 
playing a role, but also virtual black holes. Unfortunately, however, the 
theory of gravitation is unrenormalizable. Although this difficulty prevents 
one from doing reliable calculations, but it is not unreasonable to suppose 
that a semi-classical calculation will give reasonable answers. 

Einstein gravity is controlled by the action [19] 

where G is Newton’s constant and R is the Ricci scalar for the metric Qab, 
which is dehned on the spacetime M. The spacetime boundary dM has the 
induced metric hab- The quantity K is the trace of the second fundamental 
form on the boundary dM and C[h] is a functional of h dehned so that the 
action of Minkowski space vanishes. Extremization of this action for hxed 
metric on the boundary leads to the Einstein equations for Qab in M. 

The path integral for gravity is 

J (19) 


where the integral is taken over all metrics g. Our goal is to investigate how 
black holes contribute to any amplitude in quantum gravity. We hrst assume 
that this integral can be approximated by the usual Euclidean continuation 
[6, 13, 19]. The action for a single black hole of mass m is then given by 

47rm^ 


Ii = 


M/ ’ 


( 20 ) 


where Mpi is the Planck mass. Ignoring interactions between the black holes, 
we hnd the action for a collection of N black holes to be 

AttNtiR 


In — 


M/ 


( 21 ) 


In the path integral, the black holes are indistinguishable; each is indepen¬ 
dent of the others and can be positioned anywhere in space. Since N is 
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undetermined, we can evaluate Z in a box of volume V to obtain the result 

POO °° 1 T/ ^ 

dm ^ exp[-47riVm7M/] — . (22) 

The factor of V comes from accounting for the black holes being anywhere 
in the box, and the factor of 1/iV! arises from their indistinguishability. 

The combination of these results thus dehnes a probability distribution 
for having N black holes with mass m. Elementary calculations yield the 
corresponding expectation values for the number density of black holes and 
for the black hole mass, i.e., 

«pi 

where £pi is the Planck length. Thus, spacetime must be hlled with tiny 
Planck mass black holes with a density of roughly one per Planck volume. 
These microscopic virtual black holes will live roughly for one Planck time. 
This picture of the spacetime vacuum is sometimes called the spacetime foam. 

Since black holes do not conserve baryon number, these virtual black 
holes contribute to the rate of proton decay due to the gravitational interac¬ 
tion. A proton can be considered to be a hollow sphere of radius 10“^^ cm 
that contains three (valence) quarks. Suppose that two of these quarks fall 
into the same black hole at the same time. Since the black hole will decay 
predominantly into the lightest particles consistent with the conservation of 
charge and angular momentum, this process effectively converts the quarks 
into lighter particles. These particles will usually be electrons and neutrinos 
and hence baryon number conservation is generally violated. In other words, 
quantum gravity introduces an effective interaction leading to processes of 
the form 

q + q ^ i + u . (24) 

This interaction can be regarding as a four-Fermi interaction whose coupling 
strength is determined by the Planck mass. This process is mediated by black 
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holes and can violate conservation of baryon number. (Note that this process 
cannot be mediated by just gravitons because such interactions conserve both 
electric charge and baryon number.) 

The probability of two quarks being within one Planck length (£pi ~ 
10“^^ cm) of each other inside a proton is about (mp/Mpi)^ ~ 10“^^. This 
value represents the probability per proton crossing time rp ~ mj} ~ 10“^^ 
yr, if we assume that the particles move at the speed of light. In order 
for an interaction to take place (such as 0 ). a virtual black holes must 
be present at the same time that the two quarks are sufficiently near each 
other. Including this effect reduces the overall interaction probability by an 
additional factor of mp/Mpi. Converting these results into a time scale for 
proton decay, we hnd an estimated proton lifetime of rp 10^5 yr. Not 
surprisingly, this lifetime is also what one would expect from equation |]I| if 
the unihcation mass Mx is taken to be the Planck mass Mp\. 


D. Time Scales for Evaporation vs Accretion 

Once a black hole exists within a star, two very different fates are possible. 
If the black hole is large enough, it will live long enough that it accretes 
additional material from the star before evaporating away. In this case, the 
black hole can eventually accrete the entire star. On the other hand, if the 
black hole is sufficiently small, it will evaporate before interacting with the 
stellar material. In this case, the legacy of the black hole is to leave behind 
its Hawking radiation products. Since black holes are known to not conserve 
baryon number in their evaporation processes [6, 8, 13], the net result of this 
latter process is a mechanism for baryon decay. 

In order to determine the mass of a black hole required to survive, rather 
than evaporate, we set the Hawking evaporation time equal to the interaction 
time (both time scales depend on the mass Mbh of the black hole). The 
evaporation time is given by 

Tevap = T-pi(Mbh/Mpi)^ , (25) 
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where Tpi is the Planck time. The interaction time is given by 


Tint = {nav)-^ = g(Mbh/Mpi)-2(Mpi/A)3, (26) 

where we have introduced an energy scale A which dehnes the number density 
of the star through the relation n = (notice that A ~ 200 MeV for a 
neutron star and A ~ 0.1 MeV for a white dwarf). Notice also we have 
taken v = c = 1, which results in the shortest possible interaction time. The 
critical mass black hole, the smallest black hole that will survive rather than 
evaporate, is thus given by 


Mbh* = Mpi (47r)-'/5(Mpi/A)3/5 . (27) 

For nuclear densities, such as those encountered in a neutron star, the critical 
mass black hole is about Mbh* = 4 x lO^^Mpi 4 x 10® g (about the mass 
of a tractor trailer). For the density of a white dwarf, the critical mass is 
somewhat larger, Mbh* = 4 x lO^^Mpi. 

The corresponding values of N, the number of particles required to make 
the black hole, are also large: N ^ 4 x 10^® for the neutron star case and 
N ^ 4 X 10^^ for the white dwarf. 

All of the mechanisms for black hole formation favor small black holes 
over large ones. In particular, the most probable black holes will be much 
smaller than those required to survive and accrete the entire star, rather than 
evaporate. Even if the minimum possible mass for a black hole is as large as 
the Planck mass Mpi, almost all black holes will still evaporate long before 
they can interact and grow. The net result of black hole processes within 
degenerate stars is thus to create a channel for the evaporation of mass energy. 
In other words, black hole processes lead to the decay of baryons and the loss 
of energy from the star. 
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III. LONG TERM STELLAR EVOLUTION 


Given the results obtained above, the long term evolution of cold degen¬ 
erate stars is now clear: They will slowly evaporate. 

The rate at which the stars evaporate depends on the black hole formation 
rate, which in turn depends on the minimum mass black hole that can form. 
If the smallest black hole mass is less than or comparable to the proton mass, 
then the net result will be gravitationally driven proton decay at a rate Fp. 
If the minimum black hole mass is larger and contains the mass equivalent 
of N nucleons (e.g., if the Planck mass is the minimum black hole mass, N 
~ 10^®), then we can still define an effective nucleon decay rate through the 
relation 

Fp = VFbh , (28) 

where Fp is the rate of black hole production. This form does not apply 
when the black hole mass Mbh = Nmp is larger than the critical mass Mbh* 
required to survive rather than evaporate (see equation p7|). 

Once the decay rate for its constituent nucleons is determined, the star 
will evolve is a simple manner [1, 3]. The total luminosity is given by 

L* = FpM*(f) = FpM*oe-^-*, (29) 

where M* is the stellar mass and M*o is its initial mass. 

The decay products from black hole evaporation can be divided into two 
general classes: those that interact with stellar material before leaving the 
star (optically thick components) and those that leave with out interacting 
(optically thin components). Photons and all charged particles will be op¬ 
tically thick, whereas neutrinos and gravitons will be optically thin. The 
fraction T of optically thick components will generally thermalize within the 
stellar interior and ultimately produce a photon luminosity 

= ATTRlapT^ ■ (30) 
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In the second equality we have dehned the corresponding surface temperature 
of the star through the usual black body relation, where as is the Stefan- 
Boltzmann constant. The stellar radius i?* is determined by the usual mass- 
radius relation for cold degenerate stars [7, 9], 



where is the mass of the particle responsible for the degeneracy pressure, 
i.e., is the electron mass rrie for white dwarfs and the neutron mass rriN 
for neutron stars. 

Combining equations and 0. we obtain an expression for the evo¬ 
lutionary tracks of the stars in the Hertzsprung-Russell (H-R) diagram. 


T 


1.6 (.FFp)"/' ( iM±) . (32) 

^rriDmp' 

In this regime, stellar evolution is thus determined once the effective nucleon 
decay rate Tp is specihed. Unfortunately, as discussed above, the decay rate 
rp can take a rather wide range of values. The resulting evolution is shown 
in Figure 1 for the representative value of Tp = (10^® yr)~^. 

Degenerate objects become larger in radial size as they lose mass. White 
dwarfs will follow the evolutionary tracks given by equation until they 
expand to a size for which they are no longer degenerate. At this point 
in the star’s future evolution, when the stellar mass has fallen to ~ 10“^ 
Mq, the star’s track in the H-R diagram will change its slope. The above 
expressions for the stellar luminosity and the surface temperature remain 
valid, but the mass-radius relation H no longer applies. The structure 
of non-degenerate matter is largely determined by Coulomb forces, which 
enforce roughly uniform density po, and the star follows an evolutionary 
track [1] given by 

.‘tfiTTlT?, 

(33) 




367ra| 
:f^Tp^pI * 
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where po ~ 1 g/cm^. At the start of this phase of evolution, the object has 
a mass, radius, and density closely akin to those of the planet Jupiter. The 
object then continues to lose its mass until it can no longer be considered a 
star. The end point of stellar evolution occurs when the mass has fallen to 
about 10^"^ g, the point where the object can no longer completely thermalize 
its internal radiation. 


Neutron stars follow a parallel evolutionary track given by equation [p2 


where m^) is now the neutron mass. However, as they lose mass and the 
neutrons are lifted out of degeneracy, neutron stars must eventually expe¬ 
rience a signihcant readjustment when the mass drops below ~ 0.0925Mq 
[7]. During this readjustment phase, a neutron star will either expand to a 
white dwarf conhguration [1] or explode and widely disperse its remaining 
mass [20]. These possibilities are illustrated by the dashed curve and arrows 
shown in Figure 1. 

For completeness and comparison, we note that stellar black holes evolve 
in a completely different manner. As black holes lose mass through Hawking 
evaporation, they become smaller while their surface temperature and lumi¬ 
nosity increase. Black holes follow an unambiguous track in the H-R diagram 
given by 

ibh = Ant ■ (34) 


A portion of the black hole evolutionary path is shown in the lower left part 
of Figure 1. This evolution continues until the black hole mass dwindles to 
the Planck mass and the surface temperature reaches the Planck tempera¬ 
ture. Beyond this critical juncture, further black hole evolution depends on 
quantum gravity effects and is largely unknown. 

Notice that the long term evolution of white dwarfs is completely specihed 
by equations 0 and [^], and that all of the parameters are known except 
for the value of the effective proton decay rate Pp. For black holes, the 
long term evolution is given by equation ||^ with absolutely no unknown 
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parameters. 


IV. SUMMARY 

The final demise of white dwarfs and neutron stars will take place on 
a time scale that vastly exceeds the current age of the universe. However, 
in the absence of conventional unification scale proton decay, the black hole 
processes outlined here will determine the ultimate fate of these cold degen¬ 
erate stellar remnants. Within the stars, the formation of small black holes 
is strongly favored over large black holes and hence essentially all processes 
result in black holes that evaporate rather than accrete additional stellar 
material. The net result is that degenerate stars will themselves evaporate 
and follow the evolutionary tracks depicted in Figure 1. Since the vast ma¬ 
jority of stars are destined to end their conventional nuclear burning lives as 
white dwarfs, with most of the remainder becoming neutron stars, this work 
describes the long term fate and evolution of almost all of the stars in the 
sky. 
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Figure Caption 

Figure 1. The long term evolution of cold degenerate stars in the H-R dia¬ 
gram. Upon completing the early stages of stellar evolution, white dwarfs and 
neutron stars rapidly cool to an equilibrium temperature dictated by black 
hole induced proton decay (assumed here to occur at rate Pp = [10"^^yr]“^). 
The white dwarf models are plotted at successive twofold decrements in mass. 
The mean stellar density (in Log[p/g]) is indicated by the gray scale shading, 
and the sizes of the circles are proportional to stellar radius. The relative 
size of the Earth and its position on the diagram are shown for compari¬ 
son. The evaporation of a 1 Mq neutron star is illustrated by the parallel 
sequence, which shows the apparent radial sizes greatly magnihed for clarity. 
The Hawking radiation sequence for black holes is also plotted. The arrows 
indicate the direction of time evolution. 
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